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ABSTRACT 

We present a 1.3 mm continuum survey of protoplanetary disks in the 2-3 Myr old cluster, IC348, 
with the Submillimeter Array. We observed 85 young stellar objects and detected 10 with 1.3 mm 
fluxes greater than 2 mJy. The brightest source is a young embedded protostar driving a molecular 
outflow. The other 9 detections are dusty disks around optically visible stars. Our millimeter flux 
measurements translate into total disk masses ranging from 2 to 6 Jupiter masses. Each detected 
disk has strong mid-infrared emission in excess of the stellar photosphere and has Ha equivalent 
widths larger than the average in the cluster and indicative of ongoing gas accretion. The disk 
mass distribution, however, is shifted by about a factor of 20 to lower masses, compared to that in the 
~ 1 Myr old Taurus and Ophiuchus regions. These observations reveal the rapid decline in the number 
of small dust grains in disks with time, and probably their concomitant growth beyond millimeter 
sizes. Moreover, if IC348 is to form planets in the same proportion as detected in the field, these faint 
millimeter detections may represent the best candidates in the cluster to study the progression from 
planetesimals to planets. 

Subject headings: stars: formation - stars: pre-main sequence - circumstellar matter - protoplanetary 
disks 



1. INTRODUCTION 

Dusty disks persist around most low mass stars for a 
Myr or more, much longer than the lifetime of the na- 
tal molecular core. At the stage where the core is dis- 
persed, the now optically visible central star has nearly 
achieved its final mass and the disk is no longer proto- 
stellar, but it may still be protoplanetary. In the liter- 
ature of young stell ar objects (YSOs), th is corresponds 



to physical sta ge II (iRobitaille et al.ll2006[ ) and observa- 
tional Class II (|Lada 
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The presence of protoplanetary disks is most readily 
inferred via optically thick infrared emission from the 
warm dust close to the star. The capacity of a disk to 
actually form planets, however, requires knowledge of its 
surface density (or at least mass) which necessitates ob- 
servations at millimeter wavelengths where the thermal 
dust emission is optical ly thin and all radii contribute 
(|Williams fc Ciezall201in . 

Compared to infrared observations, there are relatively 
few millimeter surveys of YSOs on account of the rapid 
decline of dust emission with w avelength. The pioneering 
work of iBeckwith et al.l ()1990t ) showed that many proto- 
planetary disks in Taurus had masses comparable to or 
in excess of the minimum mass solar nebula (MMSN; 
0.01 Af© ~ lOMjup). This result was confirmed and ex- 
tended t o the similarly close Ophiuc hus star forming re- 
gion by lAndre fc Montmerld (|1994l ). Improvements in 
detector technology led to a significantly increased sen- 
sitivity and the ability to detect disks with masses < 
IMjup ([Andrews fc Willianiil2005[ ). The high resolution 
afforded by interferometry can distinguish the millime- 
ter emission from individual disks in the crowded envi- 
rons of the mor e dista. nt Ori on Trapezium Cluste r where 
iWilhams et all (pOOl and lEisner fc Carpented ([20061 ) 
found several disks with the capacity to form a solar sys- 
tem. In a more comprehensive survey, iMann fc Williams 



([20091 ) showed that the Trapezium Cluster disk mass dis- 
tribution is similar to Taurus and Ophiuchus but trun- 
cated at the high mass end, > 34Mjup, due to photoe- 
vaporation by the 40 M© 06 star, 9^ OriC. 

The Taurus, Ophiuchus, and Orion star forming en- 
vironments are all quite young, with typical YSO ages 
^ 1 Myr. A key question is how does the disk mass 
evolve with time. A h andful of o lde r, m assive disks ex- 
ist (e.g., TW Hvdrae. Wilner et al.ii2000l) but the survey 
of 3 Myr to 3 Gyr stars bv lCarpenteT "etaLl ([2001 show 
that such objects are very rare. This study of IC348 is 
motivated by the desire to quantify the decline in disk 
masses after 1 Myr. 

IC348 is located at the edge of the Perseus molecu- 
lar cloud and has been well studied at optical to in- 
frared wavelengths. It contains 283 spectroscopically 
confir med members in a re latively compact (20' x 20') 
region ([Luhman et al.lll998[ ) with' a mean age estimated 
to lie between 2 and 3 Myr ([Herbid [19981 ) ■ About haff 
of the YSOs in IC348 show mid-infrared emission in ex- 
cess of the stellar photosphere ( Lada et al . 200Q) , which 
is consistent with the median disk lifetime derived fro m 
t he statistics of ma ny clusters ([Hernandez et al.ll2007[ ). 

ICarpenteil ([2002D mapped the central region of IC348 
containing 95 YSOs at 3 mm with the OVRO interfer- 
ometer but did not find any significant detections at a 
3cr sensitivity of 25 Mjup. By selecting specific fields 
centered on concentrations of known YSOs, rather than 
a contiguous mosaic, and by observing at the signif- 
icantly shorter wavelength of 1.3 mm, we are able to 
achi eve an order of m agnitude higher mass sensitivity 
than lCarpenteii ([20021 ) and make the first mass measure- 
ments of protoplanetary disks in this region. We describe 
the observing strategy in [J2| and the results in [J3l We 
analyze the properties of the detected sources and com- 
pare the distribution of disk masses with the younger 
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Taurus, Ophiuchus, and Orion regions in gl then dis- 
cuss the imphcations of this work for disk evolution and 
planet formation in Sj5] We summarize our findings in 



2. OBSERVATIONS 

We observed 22 fields containing 85 YSOs in IC348 
with the Submillimeter Array (SMA) in its compact con- 
figuration over 11 nights from October 2009 through De- 
cember 2010. The field locations were chosen to maxi- 
mize the total number of YSOs in the survey while avoid- 
ing localized regions of cloud emission deter mined from 
850 /im SCUBA maps (|Hatchell et aUlMl . The ob- 
served fields are overlaid on a three-color infrared image 
of the region created from the Spitzer Heritage Archivq3 
in Figure [T] 

Observations were generally carried out in average 
Mauna Kea weather with precipitable water vapor levels 
ranging from 1.5 to 3 mm. IC348 transits nearly over- 
head in the sky providing long observational tracks and 
low median system temperatures, T^ys = 90 — 140 K. Typ- 
ically observations of 3 fields were interleaved with each 
other on a single track so that the total time per field was 
2-2.5 hours, resulting in an rms noise ranging from 0.4 
to 0.9 mjy with median 0.6 mJy. The field centers, dates 
of observation, and rms noise are tabulated in Table [TJ 

The receivers were tuned to place the 230.538 GHz 
(1.3 mm) J = 2 — 1 line of CO in the upper sideband 
and the same transition of the isotopologues, ^^CO, and 
C^^O in the lower sideband. The line- free regions pro- 
vided nearly 4 GHz of bandwidth for continuum measure- 
ments. 

Together with the observations of IC348, calibra- 
tion measurements were made through observations of 
3C454.3, 3C279, or 3C273 for bandpass, J0336-H323 and 
3C84 for amplitude and phase, and Uranus and Titan for 
the absolute fiux scale. The data were calibrated using 
standard procedures in the MIR software packagfl The 
rms amplitude and phase noise in the calibrated data 
were 20% and 10° - 20° respectively. The calibrated 
visibilities were then exported into MIRIAE^Ifor imaging 
and analysis. As protoplanetary disks are too small to 
be resolved in these compact configuration data, the im- 
ages were produced using natural weighting to produce 
the highest signal-to-noise and an approximately circu- 
lar synthesized beam ~ 2'.'5 beam. The flux density or 
upper limit toward each YSO was determined by fitting 
a point source with the appropriate phase offset in the 
visibility plane. 

3. RESULTS 

Contour maps of the 1.3 mm emission toward each 
of t he 22 observed fields are in the (online) Fig- 
ures ~ . T ~. 



2a [Online-only] I2bl These fields c ontain a total of 



85 Y SOs trom the iLiihman et ai.l ()2003f ) compilation, all 
of which were det ected with Spitzer in multiple bands 
and tabulated by ILada et all (|2006D . We detect 10 of 
these YSOs with a signal-to-noise ratio greater than 3. 
These are labeled in the figures with the identification 

^ http: / /archive. sp itzer. caltech.edu/ 
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Figure 1. Three-color infrared image of IC348 showing YSOs and 
nebulosity of the associated molecular cloud. The Spitzer MIPS 
24 fim observations are in red, IRAC 5.8 ^tm in green, and IRAC 
3.6 /im in blue. The blue circles show the location and 30" extent 
of the 22 SMA fields. The alphabetical labeling follows Table [T] 
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number from Table 1 of ILada et all (|2006[) for ease of 
reference to the mid-infrared fluxes later. 

There are also a small number of possibly significant 
emission peaks that do not correspond to the position 
of a known YSO that may be due to fluctuations in the 
background cloud emission or sidel qbes from strong e mis- 
sion outside of the primary beam ([Hatchell et al.| [2005'). 
Strong CO and, in many cases, isotopologue emission is 
also found but the images do not show clear evidence 
for point source emission toward the YSOs and we are 
unable to constrain the gas content of the disks. Cloud 
contamination is much more severe in the line observa- 
tions, possibly due to the emission being optically thick 
and perhaps more spatially variable than the continuum. 
A similar s ituation was found in the Trapezium Cluster 
(jMann fc W illiams 2010). 

The strongest source, by far, is a core that is appar- 
ent as an isolated p oint source in the SCUBA map of 
iHatchell et al.l (|2005[ ). This is clearly a less evolved, pro- 
tostellar source and we describe it briefly in the follow- 
ing subssection for completeness. As the emission in this 
source is dominated by the protostellar envelope, we do 
not include it in the discussion and analysis of the pro- 
toplanetary disks thereafter. 

3.1. The Class I Protostar, IRAS 03^10+3152 

This bright object was chosen to lie in the first field as a 
simple test of the observational setup. It is a well known 
Class I YSO, cata.loged as IRAS 03410 +3152, source 
51 in ILada et al.l (|2006D . source 101 in iHatcheh et"all 
((2005|1 . and source 29 in the Chandra X-ray study of the 



Protoplanetary disks in IC348 



3 
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Figure 3. The ^'^CO{2 - 1) outflow in the Class I source, IRAS 
03410-1-3152. The blue contours show the integrated intensity over 
the velocity range 4.5 — 6kms~^, while the red contours show the 
intensity integrated over 11-19 kms-l. The black cross marks the 
position of the continuum emission peak. Red contour levels start 
at 0.8 Jy and increase by multiples of 2.4 Jy, while blue contour 
levels start at 0.2 Jy and increase by multiples of 0.2 Jy. 



cluster by iPreibisch &: Zinneckeii ()2002f) . It has a bolo- 
metric temperature of 4 63 K and luminosity of 1.6^0. 
iWalawender et all (|2006[ ) also detect several H2 shocks 
on either side of the star. 

We detect both strong continuum and line emis- 
sion toward IRAS 03410+3152. The continuum is 
a p oint source at the resolution of our data (Fig- 
|2a [Onhne-onlyJI but the CO 2 
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prominent red- blue outflow lobes with a moderate open- 
ing angle (Figure El) characteristic of Class I sources 
(jArce fc Sargentl[20Cll . 

IC348 abuts an active star forming region in the 
Perseus cloud. IRAS 03410+3152 belongs to this younger 
population but the SMA field contained a second source 
(Lada source ID 468) that is a bona-fide Class II YSO. 
Our detection of this object in the same field provides 
a graphic illustration of the tremendous difference in 
the millimeter fiux from a envelope-dominated to disk- 
dominated YSO. The rest of this paper concerns the na- 
ture of the latter objects. 

3.2. Protoplanetary Disk Masses 

Nine of the remaining 84 YSOs in the 22 SMA fields 
were detected with 1.3 mm fluxes ranging from 2 to 
6.5 mJy. The 3ct limits on the fluxes of the 75 non- 
detections were typically less than 2mJy. We stacked 
the emission of these non-detections to constrain their 
average properties. There was no clear detection in 
the stacked map with a 3cr limit to the average fiux of 
0.3 mJy. 

As the millimeter emission is optically thin, it provides 
a direct measure of the dust mass. We convert fluxes to 
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where the dust opacity, K230GHz = 0.023 cm^ g~^, follows 
the commonly used prescription from ' Beckwith et al.l 
([1990) and implicitly assumes an ISM gas-to-du st ratio 
of 100 . We use a distance d — 320 pc based on iHerbi"^ 
(|1998f ). Although the dust in the disks lies over a range 
of radii and therefore temperatures, the Planck func- 
tion, B^, is calculated at a a single temperatur e 20 K 
based on the modeling by [Andrews fc Williamsl ()2005l ) 
who showed this provides an excellent fit to detailed mod- 
eling of the infrared-millimeter spectral energy distribu- 
tion (SED). 

The vernacular of disk mass, as derived in this way, is 
standard in the literature but is really a shorthand for a 
procedure that has much more validity in the interstellar, 
as opposed to circumstellar, medium. It is important to 
keep in mind that we are really only estimating the mass 
of micron to millimeter sized particles and then extrap- 
olating the gas content by multiplying by two orders of 
magnitude! 

The inferred disk masses vary between 2 and 6Afjup. 
Their properties are tabulated in Table [1] The survey is 
100% complete for disk masses Md > 2.3 Mjup. The limit 
on the stacking of the non-detections implies that the 
average mass of the non-detections, (Md) < 0.27 Mjup. 

The nine detections, although only a small fraction of 
the entire sample, provide the first disk mass measure- 
ments in IC348. They provide important constraints 
on the disk evolution at 2-3 Myr with commensurate 
implications for understanding planet formation. We 
place these measurements in context with disk masses 
in younger star forming regions below after first com- 
paring the properties of the detected and non-detected 
sources. 

4. ANALYSIS 

4.1. Comparison of Detected and Non-detected Sources 

Figure |4] plots the SED of the nine detected sources. 
The mid-infrared fiuxes in Spitzer IRAC bands 1-4 
(3.6 - 8.0 Mm) a nd MIPS band 1 (24 //m) are from 
iLada et al.l (|2OO60 . Optical and near- infrared flux den- 
sities in the i?c, ^C, J, H, and ifg-bands we re obtained 
by cross-referencing the source positions in ILada et al.l 
i200(l with tables p ublished bv Icieza et all (|2007| ) a^ 
fLuhman et al.l (|2003D . We estimate the extinction. Ay = 
A. imR - I) = 4.76[(i?c - /c)obs - [Re - IcH as 
in iCieza et al.l 1)2007 ') where the intrinsic stellar colors, 
{Rc — Ic)o, are from Kenyon & Hartmann (1995) for the 
spectral types as determined bv ILuhman et al.l (|2005| l. 
The infrared SEDs are similar, but marginally lower in 
a handful of cases, than the median SED of classical T- 
Tauri stars (CTTS) out to 24 /im, which is over-plotted 
in each panel. 

Many of the millimeter non-detections have quite sim- 
ilar infrared SEDs as the millimeter detections in Fig- 
ure m However, there is a statistical difference in 
the accretion properties of the two samples. Figure [5] 
plots the cumulative distribution of the Ha equivalent 
width s, Wa(Hq;), fo r 60 o f the 84 disks in our survey 
from iLuhman et all (|2003f ). All nine disks detected at 
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Figure 4. Spectral Energy Distributions for ttie nine SMA de- 
tected disks labeled by the IDs in Table [l] The blue squares repre- 
sent the observed optical and IRAC fluxes, uncorrected for extinc- 
tion, while the filled red circles represent the extinction-corrected 
optical and IRAC fluxes, as well as the MIPS 24 /im and SMA 
1.3 mm fluxes. Purple diamonds represent upper limits to the 
MIPS 24 ^m flux, and the black crosses represent l-cr error bars. 
The green lines show model photospheres for the published spectral 
types of each star, and the black solid and dashed lines show the 
median and quartiles, respectively, of the infrared SEDs for CTTS, 
scaled to the J-band flux. 




W^(Ha) [Angstroms] 

Figure 5. Cumulative distributions of Ha equivalent widths for 
the SMA detected (red line) and non-detected sources (blue line), 
for the 60 out of 84 surveyed disks where such measurements exist. 



1.3 mm have large Wa(Hq!), indicative of high accre- 
tion, but most of the non-detections have much lower 
values. About two-thirds of the non-detections with 
WOi(Ha) measurements have equivalent widths less than 
the minimum of the detected disks. Based on the spectral 
type d ependent W\ (Ha) cutoff defined by White & Basri 
(|2003). most of the sources that are not detected in our 
SMA survey are weak- lined T-Tauri stars (WTTS). How- 
ever, we only detect 9 out of the 25 identified CTTS. 
The sensitivity limitations of our survey prevent us 
from drawing any strong conclusions, but the general 
correspondence between the millimeter detections and 
high WA(Ha) suggests that the gas and dust depletion 
timescales are not greatly dissimilar. 

The situation is encapsulated in the IRAC color-color 
plot in Figure ini where different symbols represent CTTS, 
WTTS and SMA detections. The vertical a nd horizon- 
tal dotted lines show the boundaries used bv iCieza et al.l 



(|2007D to distinguish between YSOs with colors consis- 
tent with a stellar photosphere (lower left corner), tran- 
sition disks with excess infrared emission only apparent 
at wavelengths greater than 4.5 /im, and "full" disks with 
excesses at 4.5 and 8 fim. All the CTTSs in our sample 
display excess emission at 8 /im, and the nine disks de- 
tected by the SMA exhibit excess emission at both 4.5 /im 
and 8 /Ltm. However, most disks with an infrared excess 
were not detected in our millimeter survey and there is 
no correlation between disk mass and the amount of in- 
frared excess. We stacked the SMA emission from all of 
the stars with infrared excesses that were not individu- 
ally detected with the SMA, but did not find a significant 
detection. 

T here are 45 Chandra X -ray sources within our sur- 
vey (jPreibisch fc Zinneckerl [2002'l. 40 of these are coin- 



cident with known cluster members from iLuhman et al.l 
liOOS) and are active YSOs. 5 were detected with the 
SMA (including IRAS 03410-^3152) but we do not see 
any clear trend between the millimeter properties of the 
X-ray detected and undetected sources. The other 5 
Chandra sources have no optical or infrared counter- 
part to a level that rules out highly extincted substel- 
lar o bjects. These are in dicated by blue crosses in Fig- 
_2a [Online-onlyj|2bl and, as iPreibisch fc Zinneckeii 
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(|2002l) concluded, are most likely extragalactic contam- 
inants. Interestingly, we detect one of these, source 195 
in field V (the lowest noise map in our survey), with a 
flux 1.2±0.4mJy. 

Most YSOs in IC348 are between 2 and 3 Myr old but 
there is likely to be a significant dispersion. A simple 
explanation of our results would be that the detected 
disks are younger than the non-detections. In principle, 
protostellar ages and masses can be inferred by compar- 
ing their luminosity and effective temperature against 
model isochrones of pre-main-sequence tracks but, in 
practice, the associated errors are large. Usin g the model 
isochrones of iD'Antona fc Mazzitellil ()1997[ ). we do not 
find any significant difference in the estimated ages of 
the millimeter detected and non-detected YSOs. This is 
not too surprising given previous failed searches for dif- 
ferences between CT TS and WTTS in individual clusters 
(jComez et al.|[T99^ which show that age is not the sole 
determinant of disk evolution. 

Protoplanetary disk masses are also k nown to scale 
with stellar mass (| Williams fc Ciezall20TTI ). The inferred 
stellar masses of the detected disks are all sub-solar but 
range widely from 0.02 to 0.77 Mq. In this small sample, 
there is no evidence for a dependency on stellar mass but 
the detection of a disk around source 468, a M8.25 brown 
dwarf, is noteworthy. 

4.2. Comparison with Taurus, Ophiuchus, and Orion 

Based o n the lack of det ections of relatively mas- 
sive disks, iCarpenteil (|2002f ) concluded that disks in 
IC348 have significantly lower masses than those in 
Taurus. The high sensitivity of this SMA survey has 
allowed us to make the first mass measurements of 
IC348 disks and the results provide an important bench- 
mark for tracking disk evolution at timescales compa- 
rable to their statistical half-life. Due to the large 
number of non-detections and because the sensitivity 
of the observations varied from field to field, we plot 
the cumulative disk mass distribution, -F(< m), in Fig- 
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Figure 6. Infrared color-color plot of the surveyed YSOs from 
Spitzer IRAC bands 1, 2, and 4. The abscissa is the magnitude 
difference between 4.5 and 3.6 /im and the ordinate is the magni- 
tude difference between 8.0 and 3.6 fim. The dashed lines show the 
boundaries where there a clear infrared excess above the stellar 
photosphere. The symbols characterize the sources through their 
Hcf equivalent width and millimeter emission. The SMA detections 
are all accreting and have infrared excesses at both 4.5 and 8.0 ^m. 



ure m using: the Kaplan-Mei er product limit estima- 
tor fFeig elson &: NelsonI I1985D . Similarly defined dis- 
tributions are plotted for the (sub-)millimet er derived 
disk m ass measureme nts in Taurus ([An drews fc Williams! 
120051 ). Ophiuc hus ("Andrews fc wTlliama i2007[) . and 
Orion (jMann fc Williams 2010). 

As noted by Andre ws fc Williamj (2007^, the Taurus 
and Ophiuchus distributions are remarkably similar. Be- 
cause of the greater distance to Orion, the completeness 
limit is higher but also shows a similar distribution be- 
tween 3 — 10 Mjup \ yith a decline at higher m asses due 
to photoevaporation (|Mann fc Williamsll20l0l) . A more 
conventional, binned distribution of the Taurus, Ophi- 
uchus and IC348 disk masses is plotted in Figure [H 

The disk mass distribution in IC348 is obviously very 
different from these three younger star-forming regions. 
Our inferred masses, which are the high mass tail of the 
full distribution, are comparable to the median mass in 
the other regions. Or, to put it another way, more than 
60% of Taurus and Ophiuchus disks, but less than 10% of 
IC348 disks, have masses greater than 2Mjup. Further, 
although the number of detected disks is small, there 
is no sharp falloff at the high mass end of the IC348 
distribution as in Orion. Rather, the IC348 disk mass 
distribution seems to be systematically shifted down by 
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Figure 7. Cumulative distributions of protoplanetary disk masses 
in IC348, Taurus, Ophiuchus, and Orion. The mass distribution 
for non-accreting Class III sources in Taurus is shown as the blue 
dashed line. The distributions were calculated using a Kaplan- 
Meier product limit estimator that allows the incorporation of 3o" 
upper limits into the distributions. The long vertical lines dropping 
down to zero indicate the lowest mass or upper limit in the sample. 



about a factor of 20 relative to Taurus and Ophiuchus, 



■FlC348(< m) C± fTaurus(< 20m), 



TO > 2 Mju 



The comparison here is of Class II YSOs, or optically 
visible stars with infrared excesses and signatures of ac- 
cretion but no envelope material. A handful of Class 
III YSOs (no infrared excess or ac cretion) were detected 
in the sensitive Taurus survey by [Andrews fc Williamsl 
(|2005f ) and their distribution is plotted as the dashed 
line in Figure H The 2-3 Myr IC348 Class II disks are 
more massive than the ^ 1 Myr Taurus Class III sources. 
We conclude that we are witnessing the tail end of the 
Class II phase, just before accretion ends and the disk 
d isappears. 

iLuhman et al.l ()2003() find that the distribution of stel- 
lar masses peaks at lower masses in IC348 compared to 
Taurus. The difference between the stellar mass func- 
tions is most pronounced above 1 Mq . Almost all the 
stars used for the disk mass function comparison here 
are spectral types K and M and have sub-solar masses. 
We also compared the disk mass distribution for M stars 
only and find, with statistically less precision, the same 
large difference. The dramatically lower disk masses in 
IC348 are not due to differences in stellar mass. 

Because of the steps involved in determining masses 
from millimeter fluxes, there may be systematic ef- 
fects in the comparisons between different regions. In 
each case, however, the masses were determined assum- 
ing the same char acterist ic dust temperature of 20 K, 
the same Bcckw ith et al.l (|1990l ) dust opacity, k^, — 

0.1(j//10^^ Hz) cm^ g ^, and the same gas-to-dust ratio of 
100. The average dust temperature and frequency depen- 
dence of dust opacity can be checked with SED modeling 
and are likely to be reasonably accurate. Cloud confusion 
prevents us from making an unambiguous CO measure- 
ment of gas in the disks but using the same gas-to-dust 
ratio for all regions allows the simplest comparison. Tak- 
ing this factor out, we can state with more certainty that 
there is a reduction of about a factor of 20 in the amount 
of millimeter and smaller sized dust grains in the IC348 
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Figure 8. Mass histogram of protoplanetary disk masses in Tau- 
rus, Ophiuciius, and IC348. Tlie blue bins show data for which the 
survey in each region is complete. 

disks compared to the younger regions and, presumably, 
this is due to the additional ^1 — 2 Myr of evolution. 

5. DISCUSSION 

Sensitive infrared surveys of many clusters of differ- 
ent ages have firmly established the trend of decreasing 
disk fraction with age and a me dian disk lifetime aroun d 
low mass stars of ^ 2 — 3 Myr (|Hernandez et al]|2007[ ). 
The fraction of accreting systems declines on a similar, 
thoug h possibly slightly shorter timescale (jFedele et alj 
[2nT(T i. IC348 disks are old enough that we expect many 
to be fully depleted of both dust and gas but we have 
been able to quantify the amount of dust in a few of the 
sources with infrared excesses and found that they are all 
active accretors with Ha equivalent widths greater than 
the cluster average. 

There is a wide diversity of disk SE Ds in IC348 
(|Lada et all [20061: ICurrie fc Kemro] [200l . Based on 
Spitzer IRAC and MIPS data, they consider three cate- 
gorizations: (1) anemic disks with excess emission above 
the stellar photosphere in all IRAC bands (A = 3.6 — 
8 /im) but at levels significantly below the median disk 
SED in Taurus; (2) transition disks or disks with inner 
holes that have a mid- infrared dip in their SED, with 
strong excess emission only at A > 24 /im; (3) primordial 
or "full" disks that show excess emission at all bands at 
levels similar to the median disk SED in Taurus. At our 
survey detection threshold of about 2 mJy, correspond- 
ing to 2 Mjup under the assumptions in ii3.2[ we only 



detected a handful of primordial disks. The others, in- 
cluding most primordial disks and all anemic and tran- 
si tion disks, h ave very low masses. 

iLada et all ([2006) suggested that the low infrared ex- 
cesses in anemic disks are due to a flattened geometry 
arising from dust settling. Our observations show that 
there is a paucity of millimeter-sized dust grains, how- 
ever, and the disk has not just settled but has largely 
depleted. Similarly the low masses of the transition disks 
reveals that, in general, the inner disk only dissipates af- 
ter the outer disk is largely depleted. Cicza et al. (20081) 
came to the same conclusion based on their low limits to 
the disk masses of WTTS. This is consistent with photo- 
evaporative models of disk evolution whereby the outer 
disk feeds the inner disk until the viscous accretion rate 
drops below the photoevaporation rate from the central 
star. Our low disk mass upper limits suggest that the 
photoev aporation rate is low, c onsistent with FUV ra- 
diation (jWilliams fc Ciez"all201lL and references therein). 
In fact, the disks that we have detected are the least 
evolved, in the sense that they exhibit strong infrared 
excesses at all wavelengths and Ha emission indicative 
of gas accretion. Yet we only detect about one third of 
such CTTS and their masses are all substantially below 
the canonical 10 Mjup MMSN. 

Assuming that the giant planet frequency in IC348 is to 
match that measured in the field, the absence of MMSN 
disks implies that planet formation is either complete 
or has progressed substantially. However, the possibility 
that giant planets have already for med is hard to rec- 
oncile with core acc retion models (jPollack et al.l 119961 : 
iHubickvi et al.ll2005[ ). The formation of Jupiter and, es- 
pecially, Saturn (0.3 Mjup) mass planets requires several 
Myr, significantly longer than the median YSO age in 
IC348. 

On the one hand, our results may be considered a point 
in favor of very rapid giant plane t formation through 
gravitational instability ()Bosslll997f) . On the other hand, 
if planets form over several Myr through core accretion, 
this process must have progressed to the point where 
most of the solid material is i n large bodies and be yond 
our ability to detect directly (jOreaves Ricel[2?)Tnl) . At 
our observing wavelength, A = 1.3 mm, the survey is 
effectively bli nd to particles with sizes a > 3A ~ 4 mm 
(|Drainell2006l ). For a constant dust mass, and a grain size 
distribution, n{a) a~^'^, the millimeter emission drops 
by our measured factor of 20 once the maximum grain 
size reach es about a met e r, as verified in the detailed 
models of iD'Alessio et al.l (|2001| ). 

We suggest that the primordial disks with millimeter 
emission in IC348 may be excellent candidates for study- 
ing the progression from planetesimals to planets. These 
are the most massive disks in the young cluster but, given 
their low millimeter luminosities, most of the dust mass 
probably resides in macroscopic objects. Their contin- 
uous infrared SEDs show that no planets have dynami- 
cally cleared out radial gaps or holes that are optically 
thin to starlight and their large Ha luminosities indi- 
cate that they retain some gas. Recent calculations by 
iJohansen fc Lacerdal (|2010[) and lOrmel fc Klah^ ()2010f ) 
show that centimeter-sized "pebbles" in gas rich disks 
experience strong drag forces and rapidly accrete onto 
protoplanets thereby accelerating planetary growth. Our 
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conclusion that most of the solids in the disks have aggre- 
gated beyond millimeter sizes within 2-3 Myr is also con- 
sistent with cosmochemical measurements of the absolute 
ages of chondrules within rocky meteorites f|Amelin et alJ 

[2nnl . 

6. SUMMARY 

We have conducted a 1.3 mm SMA survey of 85 YSOs 
in the nearby cluster IC348. With a mean age of 2-3 Myr, 
the cluster presents an opportunity to study YSOs at a 
critical time in their evolution when they begin to lose a 
large fraction of their disks. Our main findings are: 

• Nine protoplanetary disks are detected with in- 
ferred masses ranging from 2 — 6Mjup. 

• The millimeter detected disks are all actively ac- 
creting and exhibit excess photospheric emission 
above the photosphere at A > 4.5 /im with infrared 
SEDs similar to the median in Taurus. 

• The protoplanetary disk mass distribution in IC348 
is shifted to lower masses by about a factor of 20 
compared to the ^ 1 Myr old Taurus and Ophi- 
uchus star-forming regions. We suggest that this is 
due to substantial grain growth such that most of 
the solid mass resides in large particles, ^ 4 mm 
and perhaps up to meters, in at least some of the 
disks. 

• Grain growth beyond millimeter sizes in gaseous 
disks may promote rapid planetary growth and the 
disks that we have detected with the SMA likely 
signpost the best candidates in the cluster for wit- 
nessing the birth pangs of giant planets. 

Future directions include resolved observations to 
see whether any protoplanets might have opened up 
"dusty gaps" that are infrared bright but millimeter-faint 
(| Andrews et al.l 1201 If) . It will also be very interesting 
to try additional spectroscopic observations to measure 
their gas content, and to observe them at longer wave- 
lengths to constrain the population of centimeter and 
larger sized grains. Due to their low millimeter luminosi- 
ties, however, these are tasks for upcoming facilities. 

We thank Anders Johansen, Jes Jorgensen, and Sean 
Andrews for interesting discussions on the topics of grain 



growth and disk evolution. This work is supported by the 
NSF through grant AST08-08144. 
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Table 1 

Observed fields 



Field 


R.A. 




Decl 


Date 


RMS 




(J2000) 


(J2000) 


fvvmmdd) 




A 


03:44:12.80 


OZ. 




091013 


0.48 


B 


03:44:44.57 




■1 fi-AQ fin 


091013 


0.51 


C 


03:44:23.45 


OZ. 




091013 


0.46 


D 


03:44:38.19 


OZ. 


.Uo.io.oy 


091020 


0.80 


E 


03:44:35.91 


OZ. 


no-fi/i no 


091020 


0.75 


F 


03:44:37.26 


OZ. 


.uo.zi.oy 


091020 


0.77 


G 


03:44:30.31 


OZ. 


no -/I 7 y1 o 

,uy.4 1 .4y 


091204, 091206 


0.73 


H 


03:44:29.52 


OZ. 


ni -1 1 oi 


091204, 091206 


0.68 


I 


03:44:40.71 


32: 


:09:56.90 


091204, 091206 


0.69 


J 


03:44:32.29 


32: 


:05:41.09 


091207 


0.62 


K 


03:44:13.64 


32: 


:13:27.60 


091207 


0.59 


L 


03:44:42.96 


32: 


:08:44.29 


091207 


0.61 


M 


03:44:37.60 


32: 


:11:57.59 


091222 


0.50 


N 


03:44:21.55 


32: 


:12:07.69 


091222 


0.49 


O 


03:44:26.39 


32: 


:08:12.49 


091222 


0.50 


P 


03:44:18.78 


32: 


:07:32.50 


091227 


0.46 


Q 


03:44:34.38 


32: 


:12:58.39 


091227 


0.48 


R 


03:44:45.34 


32: 


:04:01.79 


091227 


0.46 


S 


03:44:32.04 


32: 


:11:43.70 


101006 


0.94 


T 


03:44:43.53 


32: 


:07:42.74 


101019 


0.82 


U 


03:44:36.96 


32: 


:06:45.20 


101122 


0.91 


V 


03:44:56.15 


32: 


:09:15.19 


101214 


0.37 



Table 2 

Detected Disk Properties 



ID 


R.A. 


Decl 


Field 


Spectral Type 


Fl.3 






WxiHa) 


Mstar 




(J2000) 


(J2000) 






(mJy) 


(mJy) 


(Mjup) 


(A) 


(M@) 


(1) 


(1) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


51 


03 44 12.97 


32 01 35.40 


A 




75.90 


0.55 


69.19 


45 




153 


03 44 42.77 


32 08 33.90 


L 


M4.75 


6.49 


0.71 


5.92 


40 


0.20 


32 


03 44 37.89 


32 08 04.20 


D 


K7 


5.47 


0.81 


4.99 


68 


0.77 


221 


03 44 40.26 


32 09 33.20 


I 


M4.5 


3.23 


0.89 


2.94 


40 


0.20 


248 


03 44 35.95 


32 09 24.30 


E 


M5.25 


2.81 


0.81 


2.56 


30 


0.15 


468 


03 44 11.07 


32 01 43.70 


A 


M8.25 


2.37 


0.56 


2.16 


400 


0.02 


8078 


03 44 26.69 


32 08 20.30 


O 


M0.5 


2.07 


0.51 


1.89 


75 


0.65 


100 


03 44 22.32 


32 12 00.80 


N 


Ml 


2.06 


0.51 


1.88 


90 


0.60 


192 


03 44 23.64 


32 01 52.70 


C 


M4.5 


2.06 


0.57 


1.88 


40 


0.23 


15 


03 44 44.72 


32 04 02.70 


R 


M0.5 


1.98 


0.46 


1.80 


36 


0.65 



Note. — (1) Source ID and coordinates as given in ILada et al.1 ll200^). (2) Observ ed field, as labeled 
in Figure [T] (3) Spectral type of object where available, from iLuhman et al.l II2003I) . (4) 1.3 mm flux 
density. (5) Icr errors in fiux density. (6) Protoplanetary disk mass in units of Mj^p = 1/1O47M0. (7) 
Ha equivalent vfidths, from [Euhman et al.l II2003I V (8) Stellar Masses estimated from pre-main-sequence 
evolutionary tracks. 
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Figure 2a [Online-only] . Cleaned SMA maps of the 22 fields in the survey. The letter at the top left of each panel corresponds to the 
labels in Figure[l] The rms noise, a (mjy), in each image is given in the bottom left of each panel. Contour levels begi n at Sct and increase in 
steps of 2(T. Dotted contours show negative emission on the same scale. The locations of the cluster members from the lLuhman et aP 112003 ) 
compilation are shown by red crosses. All these sources were also detected with Spitzer by Lada et al. (2006). There are an additional 5 
sources with Chandra X-ray detections in the Preibisch fc Zinnecker, 1,2002 ) compilation. These are most likely extragalactic sources and 
one, PZ195 in field V, shows significant 1.3mm emission. The SMA detected YSOs are labeled by their ID number in Table[T] 
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Figure 2b. Fig 2 (contd.) 



